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O ■ Abstract 
(D ; 

pL| i We compute of the J/ip production in pp collisions at RHIC at ^/s = 500 and 200 GeV 

by using both the color octet and singlet models in the framework of non-relativistic QCD. The 
J/ip we compute here includes the direct J/^ from the partonic fusion processes and the J/V' 
lO . coming from the radiative decays of xj's both in the color octet and singlet channel. The high 

I Pt J/^P production cross section is computed within the PHENIX detector acceptance ranges: 



—0.35 < r] < 0.35 and 1.2 < |r/| < 2.4, the central electron and forward muon arms. 
It is found that the color octet contribution to J/V' production is dominant at RHIC energy in 



Q-ii comparison to the color singlet contributions. We compare our results with the recent preliminary 

data obtained by PHENIX detector for the high pj- J /ip measurements. While the color singlet 



model fails to explain the data completely the color octet model is in agreement with the single 



^ . data point above 2 GeV transverse momentum. A measurement of J/'0 production at RHIC in 

Cd \ the next run with better statistics will allow us to determine the validity of the color octet model 

of J/ip production at RHIC energies. This is very important because it is necessary to know the 
exact mechanism for J /ip production in pp collisions at RHIC if one is to make predictions of J /ip 
suppression as a signature of quark-gluon plasma. These mechanisms also play an important role 
in determining the polarized spin structure function of the proton at RHIC. 
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Recently there has been much progress in understanding of heavy quarkonium produc- 
tion mechanism in high energy hadronic colhsions. The two prominent models for heavy 



quarkonium production mechanism which are quite successfu 



at different collider energies 



are 1) color singlet model and 2) color octet model 0,13,14,01. The J/ip production 

singlet model for 

high pt quarkonium production |3| in that center of mass energy regime. The color singlet 
model gives a reasonable description of large pt J/ '4' production at ISR energies ^. In 
this context large pt J/ production data from pp collisions at ^/s = 500 and 200 GeV at 
RHIC might tell us the relative importance and validity of these two models in this new 
energy range. In particular a valid J/ip production mechanism at pp collisions will play an 
important role in determining the mechanism for J/ ip suppression at Au-Au collisions which 
is suggested to be a prominent signature for quark-gluon plasma detection at RHIC 0, 1^ . 
These mechanisms also play an important role in determining the polarized spin structure 
function of the proton. 

In the color singlet model the quarkonium is formed as a non relativistic bound state 
of heavy quark-antiquark pair via static gluon exchange. In this model it is assumed that 
the QQ is produced in the color singlet state at the production point with appropriate spin 
(S) and orbital angular momentum quantum number (L) to evolve to a bound state '^^'^^Lj 
with total angular momentum quantum number (J). The relative momentum of the QQ 
pair inside the quarkonium is assumed to be small compared to the mass Mq of the heavy 
quark so that the Q and Q will not fly apart to form heavy mesons. The non-relativistic 
non-perturbative wave functions and its derivatives appearing in the color singlet model 
calculation are either determined from the potential model or taken from experiment. 

On the other hand in the color octet model relativistic effects are taken into account 
which are neglected in the color singlet model. In the color octet model using an effective 
field theory called non-relativistic QCD (NRQCD), the dynamical gluon enters into Fock 
state decompositions of the quarkonium states. In NRQCD the expansion is carried out in 
terms of the relative velocity v (f^ ~ 0.23 for CC system and 0.1 for BB system) of the QQ 
bound state. The NRQCD Lagrangian density is given by: 

^NRQCD flight ~l~ ^heavy ~l~ ^correction (l) 
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where 

l^i^^^i = F;,, + <ll,D^[A]q, (2) 

is the usual QCD Lagrangian for gluons and the hght flavors and 

Ueavy = ^H^^t + + + (3) 

is the leading order heavy quark part and 
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[ip^D -E-E- D)ip + (j)\D -E-E- 



8M2 

rw.tn- . (' n V R — R ^ n\oi, m /At 



8M2 



[V^V -{Dx E-Ex D)i) + (j)'^a-{Dx E-Ex D)(f)] 

r 



+ + (4) 

is the correction term. In the above equations ■?/', (/) are the two component Dirac spinors 
of the heavy quark. The various correction terms in the above equation come from the 
corrections to kinetic energy, electric field coupling, magnetic field coupling and from double 
electric scattering etc.. In NRQCD the dynamical gluon enter into the Fock state decom- 
positions of different physical states. The wave functions of S-wave orthoquarkonia state 
\'^Q > appear as: 

>= 0(l)|ggp5«] > +Oiv)\QQ['PP]g > +Oiv')\QQ['S['%g > + 

Oiv')\QQ['sj,% > +Oiv')\QQ['D^j%g > + (5) 

and the wave functions of P-wave orthoquarkonium state \xqj >■ 

\xQj >= o(i)|ggppP] > +oiv)\QQ['s[% > +.... (6) 

In the above equation (1,8) refers to singlet and octet state of the QQ pair. After a QQ 
is formed in its color octet state it may absorb a soft gluon and transform into \xqj > via 
eq. © and then become J/ip by photon decay or the QQ in the color octet state can emit 
two long wavelength gluons and become a J/ip via eq. © and so on. All these low energy 
interactions are negligible and the non-perturbative matrix elements can be fitted from the 
experiments or can be determined from the lattice calculations. 

Using both mechanisms for heavy quarkonium production we compute the Pt distribution 
of J/ip production at RHIC energy in this paper. The J/ip production includes both direct 
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FIG. 1: Direct J /if) production at RHIC in the J/^ rapidity range —0.5 < y < 0.5 

J/ip and the J/ip coming from the radiative decay of Xj's: Xj ~^ J/ip + 'y. We will calculate 
the J/ip production in this paper by considering all the parton fusion processes. In the 
color octet model the non-perturbative matrix elements in the parton fusion processes are 
extracted in [4] by fitting CDF experimental data. 

The PHENIX detector at RHIC has measured a few data points in the pt distribution 
of J/ip production at y/s = 200 GeV pp collisions. As these data are preliminary our 
comparison to this experimental data is not very conclusive. However, in the future the 
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J/\|/ from Decays at RHIC 
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FIG. 2: J/il) production cross section at RHIC from the Xcj's decay in the J/V' rapidity range 
-0.5 < y < 0.5 



PHENIX detector at RHIC will measure high pt J j^) production up to px ~ 10 GeV 
with better statistics. The acceptance range of single electron for the PHENIX detector is 
—0.35 < rj < 0.35 and for single muon it is 1.2 < \rj\ < 2.4. We compute the transverse 
momentum distribution of J/ip production by using both the color octet and the singlet 
model in the PHENIX detector acceptance range. 

The 2 — > 2 parton fusion processes which contribute to the pt distribution of the J/ip and 
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Xcj production are order a^. The differential cross section for the production of a quarkonium 
state (^^^^Lj) with different quantum numbers (L,S,J) in the color singlet model is given 
by Q: 

X- ^{iJ - C'^'Lj)k) 5{s + i + u- M'). (7) 

TT dt 

Here Xi and xj are the momentum fractions carried by partons i and j inside the hadrons 
A and B respectively. The fi/Ai^i) and fj/B^Xj) are the parton distribution function inside 
the hadron. ^{ij — i^^^^Lj)k) is the partonic level differential cross section to form a 
quarkonium state (^^^^Lj) which involves the radial wave function and its derivatives at the 
origin. Using the delta function integration in Eq. ^ one obtains: 



da 

dpT 



{AB^J/ij,Xj + X) 



J2 J dy J dxi Xifi/Aixi,Q'^) x'jfj/Bix'j,Q 

2PT 



Xi - dt 



da , . . 



Ljk), 



where 



1 Xjy/sMre-y - 



v^i Xj^/s-MTey 
In the above M is the mass of the bound state quarkonium and 



(8) 



(9) 



/p| + M2. The 

partonic level differential cross section ^{ij — > (^^'^^Lj)k) for the gluon fusion processes are 

the quark, gluon processes are 



obtained in ll| by using helicity decomposition method anc 
derived in The {^^+^Lj)k) as given in and contain the non-relativistic 

wave function |i?(0)p (for direct J/ip process) and its derivatives |-R'(0)p (for xj processes) 



at the origin. For the non relativistic wave functions at the origin we take the Buchmu 
Type wave function with charm quark mass Mc=1.48 GeV. The numerical value is 



ler- 



|-R(0)p=0.81 GeV'^. For the derivative of the radial wave function at origin we use jli 



9 \R'm\' 



=15 MeV. 



2tt Af4 

In the color octet model the differential cross section for heavy quarkonium production 
at order is given by: 

da 



dp-. 



■{AB ijQixj) +X) = Yl J J Xifi/A{xi,Q'^) x'jfj/B{x'j,Q'^) 



Direct J/\|/ Production at RHIC 

1 < y < 2, pp Collisions 



10^ 



10' 



1/2 



s =500 GeV, Octet 



s'^'=200 GeV, Octet 



s'^'=500 GeV, Singlet 



s'^'=200 GeV, Singlet 



> 

O 10 
ft 10 

ID 
■D 

10-' 



10 



10" 



5 6 7 
Pt (GeV) 



8 



10 



FIG. 3: Direct J/ip production at RHIC in the J/ip rapidity range 1 < y < 2 



where 



2p- 



da 



ey dt 



(10) 



da 



i:j)|o> 



(11) 



In this paper we include the contribution from the ^-Pj^'* and -^Sq^^ processes which are 
present in the Fock state decomposition amplitude (0) in addition to ^Si'^ state of the color 



octet amplitude. Similar process is followed for the xj production by using © which takes 
into account the color octet '^sf"* state. These xj formed in the color octet state again decay 
radiatively to give J/ip. The total contribution to J/ip production in the color octet model 
is the sum of the direct J /if) and the contribution coming from xj's decay in the color octet 
model. In the above equation the short distance part of the partonic matrix element square in 
different quantum states for various partonic processes: T,\A{ij — * CC^^^^ L^j^]k) short\^ are 
very lengthy and are given in jj| . We use all the partonic processes in our calculation. We use 
the non-perturbative matrix elements: < Q\0^'^^^''^ (^^^^Lj)\Q > from the tabulation given 
in which are extracted from the CDF data. In the extraction of the these matrix elements 
< 0\Op(^Pii)\Q > = < 0\Ot'\^So)\0 > was set and then the values of the linear 

combination was extracted which is given by: ^"^^'^ ^2^°^^'^^ + ^'^^^^ ^ So)\o> _ (^2.2±0.5)x 
10~^ GeV^. For other matrix elements the values are given by: < 0|C8/'^(35i)|0 > = (6.6± 
2.1)10-=^ GeV^, < 0|Cf (3Si)|0 > = (9.8 ± 1.3)10-3 GeV^. We use the central values of 
these matrix elements in our calculation. 

We present our results for the rapidity range covered by the central electron arm and 
forward muon arm at the PHENIX detector at RHIC. The electron pseudo rapidity range: 
—0.35 < 1] < 0.35 corresponds to the J/ip rapidity range: —0.5 < y < 0.5. For the 
muon pseudo rapidity range: 1.2 < |?7| < 2.4 we present our results in the J/ip rapidity 
range 1 < y < 2 which is well covered by the RHIC/PHENIX experiment. So the 
results we will present in this paper is in the J/ip rapidity range: —0.5 < y < 0.5 and 
1 < y < 2. We will use GRV98 NLO Q parton distribution function throughout our 
calculation. The Q scale in the structure function and the renormalization scale are chosen to 
be Q = Mt = \l M"^ + pj., where M and pt are the mass and the transverse momentum of 
the heavy quarkonium respectively. In extracting the NRQCD matrix elements, the MRSDO 
parton distribution functions inside a nucleon was used in the color octet model calculation 
in 0|. In this paper we are using, instead the GRV98 NLO parton distribution function for 
the parton distribution functions. We have checked that there is no significant difference in 
the results obtained by using MRSDO and GRV98 in quantities calculated here. 

In Fig. 1 we present the direct J/\E' production cross section both from the color singlet 
and color octet models at \/s= 500 and 200 GeV pp collisions in the J/ip rapidity range 
—0.5 < y < 0.5. The solid line is the color octet contribution and dotted line is the color 
singlet contribution at A/i=500 GeV. The dashed line is the color octet contribution and 
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dot-dashed line is the color singlet contribution at v^=200 GeV. The direct J /ip production 
from color octet model at both the colliding energies are much higher than the color singlet 
counterparts. In what follows, we present the results of our computation for the pt range 
from 2 to 10 GeV. At smaller »t, the color singlet and octet cross sections are corrupted 

n n 

by collinear divergences as discussed in a Since the intrinsic motion of the incident 
partons inside colliding protons renders the differential cross section uncertain for < 2 
GeV, as discussed in js], ^ we have presented our results for Pt > 2 GeV in this paper. For 
low Pt calculation one might follow the procedure where one might have to take intrinsic 
motion of the incident parton into account. This calculation is, however, beyond the scope 
of the present paper. In this paper we report the predictions for high pt J / ip production. 
High Pt limit ~ 10 GeV corresponds to the possible pt measurement range at PHENIX. 
In the very high pt limit (> 6 GeV) the B decay contribution B —>■ J/ip + X and the 
gluon and charm quark fragmentation contribution might become important. However, we 
have not included these process here. In this paper we have computed only the prompt J/ ijj 
production from the parton fusion processes. These are the same contributions that were 
used in at Tevatron energies to extract NRQCD matrix elements. 

In Fig. 2 we present the predictions for the J/ip production cross section coming from 
the radiative decay of the xj's both from the color octet and singlet model at y/s= 500 
and 200 GeV pp collisions in the J/tp rapidity range —0.5 < y < 0.5. The solid line 
is the color octet contribution and dashed line is the color singlet contribution at -^5=500 
GeV. The dotted line is the color octet contribution and dot-dashed line is the color singlet 
contribution at ^/s=200 GeV. The Xcj contribution in each case in the figure is the sum of 
the contributions from Xo, Xi ^i-nd X2- The Xcj's decay coming from the color singlet channel 
is slightly higher at small pt but very much lower at large pt in comparison to the color octet 
contribution. It can also be noted from Fig. 1 and Fig. 2 that the direct J/ip production 
in the color octet channel is larger than the contributions coming from all other processes. 
Thus at RHIC energies, we find that the dominant contribution to J/ip production comes 
from this direct color octet process. 

In Fig. 3 we present the results of our computation for the direct J/tp production for both 
the color singlet and octet models in the J/-?/' rapidity range 1 < y < 2 at -yi=500 and 200 
GeV. The solid line is the contribution from the color octet process and dotted line is the 
contribution from the color singlet process at ^/s=500 GeV. The dashed line and dot-dashed 
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J/\|/ from X , Decays at RHIC 



1 < y < 2, pp Collisions 
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FIG. 4: J/il) production cross section at RHIC from the Xcj's decay in the J/V' rapidity range 
l<y<2 

lines are the similar curves at v^=200 GeV. In this rapidity regime we also find that the 
color octet contribution is much larger than the color singlet contribution. In Fig. 4 we 
present the J/ if) production cross section coming from the xj's decay in the J/ if) rapidity 
range 1 < y < 2 at ^/s= 500 and 200 GeV pp collisions at RHIC. The sohd line is the 

contribution from the color octet process and dashed line is the contribution from the color 
singlet process at v^=500 GeV. The dotted line and dot-dashed lines are the similar curves 
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at v^=200 GeV. The total XcJ contributions in the figure includes the contribution from xq-, 
Xi and X2- As ip' contribution to J/ip is small at this energy range we have not considered 
the ijj' contribution in our calculation. As can be seen again the dominant contribution to 
the J/ ip production comes mainly from the direct J/ ip production in the color octet channel. 

The total J/ip Pt distribution coming from combined color singlet model and color 
octet model contributions in all channels is plotted in Fig. 5 in the J/ip rapidity range 
—0.5 < y < 0.5 at y^=500 and 200 GeV pp collisions. The solid line is the pr distribu- 
tion of the total (octet (direct J/ip + Xj's decay) plus singlet (direct J/ip + Xj's decay)) 
J/t/j production cross section at ^/s=500 GeV. The dashed fine is the contribution from the 
color octet model and dot-dashed line is the contribution from the singlet model at the same 
energy. As pointed out in the above the total cross section is the sum of the cross section 
from direct J/ip and from the radiative decays of the xj's both in color octet and singlet 
channels. It can be seen that the total cross section is almost same as the cross section 
obtained from the color octet channels. The color singlet channel contribution is very small 
as can be seen from the dot-dashed line. 

Similar results are obtained for ^/s=200 GeV pp collisions in the same figure in the J/ijj 
rapidity range —0.5 < y < 0.5. The dotted line is the total (octet (direct J/iJj + xj's 
decay) plus singlet (direct J/ip + Xj's decay)) J/ip production cross section at v^=200 GeV. 
Long-dashed line is the octet contribution and thin-solid line is the singlet contribution. It 
can be seen from the figure that the total J/ip production cross section is almost same as 
that of the color octet model contribution and the singlet contribution is much smaller. A 
measurement of the J/^jJ differential cross section at RHIC in the future (both at ^/s = 500 
and 200 GeV pp collisions) will be able to tell us whether the color octet model reproduces 
the RHIC data. 

As PHENIX has measured few data points for pr distribution of J/ ip in the rapidity 
range in pp coUisions at ^/s— 200 GeV, we will compare our results with these available 
experimental data. As stated earlier we will compare the available data in the J/ ip rapidity 
range 1 < y < 2 as this range is well covered by the RHIC/PHENIX experiment in the muon 
channel. In Fig. 6 we present the px distribution of total J/ip production cross section from 
color singlet and octet model in the J/ip rapidity range 1 < y < 2 at y^=500 GeV. The 
solid line is the pt distribution of the total (octet (direct J/ip + xj's decay) plus singlet 
(direct J/ip + Xj's decay)) J/ip production cross section at ^/s—500 GeV. The dashed fine 
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J/\|/ Production at RHIC 
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FIG. 5: J/il) production at RHIC from parton fusion processes in the J/tp rapidity range 
-0.5 < y < 0.5 



is the contribution from the color octet model and dot-dashed line is the contribution from 
the singlet model at the same energy. The pt distribution of total J/V' production cross 
section from color octet and singlet model in the J/%1) rapidity range 1 < y < 2 at 
i/i=200 GeV pp coUisions is plotted in Fig. 6. The dotted hne is the total (octet (direct 

J/ip + xj's decay) plus singlet (direct J/ip + xj's decay)) J/ip production cross section at 
i/s=200 GeV. Long-dashed line is the octet contribution and thin-solid line is the singlet 
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J/\|/ Production at RHIC 
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FIG. 6: J /tp production at RHIC from the parton fusion processes in the J /tp rapidity range 
1< y <2 



contribution. 

Let us compare our results with the prehminary experimental data obtained at the 
PHENIX detector. As mentioned earlier, at small px the color singlet and octet cross 
sections are corrupted by coUinear divergences as discussed in jj] . Since the intrinsic motion 
of the incident partons inside colliding protons renders the differential cross section uncer- 
tain for ^ 2 GeV, as discussed in p, ^ we will compare our results at > 2 GeV with 
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the PHENIX data [15] in this paper. We compare our calculation with the only available 
data point (> 2 GeV) at pr ~ 2.5 GeV (see [l^)- The data point is presented as "+" 
in Fig. 6. It can be seen that the present experimental data fits the color octet model 
predictions of J/ip production and the color singlet model prediction is well below the data 
point. Hence the color singlet model is not the relevant mechanism at RHIC energy. As for 
the color octet model, the single data point matches with the color octet contribution, but 
many more data with high precision measurements are necessary to confirm the validity of 
the color octet model at RHIC energy. The next run results from PHENIX experiment will 
be able to shed light on this. 

In summary, we have computed the J/ip production differential cross section in pp colli- 
sions at RHIC at ^/s = 500 and 200 GeV by using both the color octet and singlet models in 
the framework of non-relativistic QCD. The J/ip production cross section we compute here 
includes the direct J/ip production from the partonic fusion processes and the J/ip coming 
from the radiative decays of xj's both in the color octet and singlet channel. The high pt 
J/'?/' production cross section is computed within the J/'?/' rapidity range: —0.5 < y < 0.5 
and 1 < |?7| < 2, in the the central electron and forward muon arms respectively. It is 
found that the color octet contribution to J ftp production is dominant at RHIC energy in 
comparison to the color singlet contributions. We have compared our results with the recent 
data obtained by PHENIX detector for Pt > 2 GeV at ^/s =200 GeV pp collisions. While 
the color singlet model fails to explain the data completely the color octet model seems to 
reproduce the data. As the present high pt data of J/ ip production at PHENIX detector at 
RHIC is very restricted, measurement of high pt J / ip production at RHIC in the next run 
with better statistics will be able to tell us the degree of validity of the color octet model 
of J/tp production at RHIC energies. It is necessary to know the exact mechanism of J/ip 
production in pp collisions at RHIC in order to make a prediction of J/ip suppression as a 
signature of the quark-gluon plasma in Au-Au collisions. Knowledge of the J/ ip production 
mechanism is also crucial for isolating the polarized spin structure functions of the proton 
from the data. 
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